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Density Functional Resonance Theory (DFRT) is a complex-scaled version of ground-state Den- 
sity Functional Theory (DFT) that allows one to calculate the resonance energies and lifetimes of 
metastable anions. In this formalism, the exact energy and lifetime of the lowest-energy resonance of 
unbound systems is encoded into a complex "density" that can be obtained via complex-coordinate 
scaling. This complex density is used as the primary variable in a DFRT calculation just as the 
ground-state density would be used as the primary variable in DFT. As in DFT, there exists a 
mapping of the ./V-electron interacting system to a Kohn-Sham system of TV non-interacting par- 
ticles in DFRT. This mapping facilitates self consistent calculations with an initial guess for the 
complex density, as illustrated with an exactly-solvable model system. Whereas DFRT yields in 
principle the exact resonance energy and lifetime of the interacting system, we find that neglect- 
ing the complex-correlation contribution leads to errors of similar magnitude to those of standard 
scattering close-coupling calculations under the bound-state approximation. 



Density Functional Theory (DFT) [IH3] provides one 
of the most accurate and reliable methods to calculate 
the ground-state electronic properties of molecules, clus- 
ters, and materials from first principles. It is one of 
the workhorses of computational quantum chemistry 0] . 
In addition, DFT's time-dependent extension (TDDFT) 
[5] can now be applied to a wealth of excited-state and 
time-dependent properties in both linear and non-linear 
regimes [5J. When the TV-electron system of interest 
has no bound ground state, however, neither DFT nor 
TDDFT can be applied in a straightforward way. A cor- 
rect DFT calculation converges to the true ground state 
by ionizing the system, thus leaving no reliable starting 
point for a subsequent TDDFT calculation on the N- 
electron system. In practice, a finite simulation box or 
basis set can make the system artificially bound [5] , 
but information about the relevant lifetimes is lost in the 
process. 

We address here this fundamental limitation of ground- 
state DFT, and propose a solution. 

Consider a system of N interacting electrons in an 
external potential v(r), with ground-state density n(r). 
The potential is set to be everywhere positive and go to 
a positive constant C as |r| — > oo. The ground-state en- 
ergy is E > 0. We start by asking how the gound state 
density changes when a smooth step is added to v(r) at 
a radius |R| that is larger than the range of v(r). The 
step is such that the new potential v(r) coincides with 
v(r) for |r| < |R| but goes to zero at infinity. Since v(r) 
is everywhere positive, all N electrons tunnel out and 
v(r) supports no bound states. The correct ground state 
energy is now E = 0, and the new density n(r) is delocal- 
ized through all space. In practical calculations, however, 
v(t) and v(r) cannot be distinguished if |R| is beyond 
the size of the simulation box. The result provided by 



ground-state DFT using the exact exchange-correlation 
functional is not E, but E > 0, and the density obtained 
is n(r) as if the system were bound. Even when the sim- 
ulation box is large enough to include the steps, use of a 
finite basis-set of localized functions will artificially bind 
all electrons. Clearly, such calculations do not provide 
approximations to the true ground-state energy and den- 
sity of u(r), but to those of its lowest-energy resonance 
(LER). 

The purpose of this letter is to establish an analog 
of KS-DFT that provides the in-principle exact LER- 
density along with its energy and lifetime for any finite 
|R|. As |R| — > oo, the results coincide with those of stan- 
dard KS-DFT. For higher-energy resonances, TDDFT is 
needed as a matter of principle [3 HH] ■ 

First, we note that as |R| — > oo, the complex density 
ng(r) associated with the LER of 

H v =f + V ee + J dvn{v)v{r) , (1) 

becomes equal to the complex density ng(r) associated 
to v(re ie ). In Eq. [l] f = -|EiIi v ? is the kinetic 
energy operator, V ee = Yrt,j& l r * _ r J'l _1 is the electron- 
electron interaction, and n(r) = YliLi ^( r — *i) 1S the 
density operator. (Atomic units are used throughout). 
To find ng(r), we complex-scale H v by multiplying all 
electron coordinates by the phase factor e tB ', diagonalize 
the resulting non-hermitian operator H®, and calculate 
the bi-expectation value of n(r) as: 

n e (r) = ^|n(r)|^) , (2) 

where l^^) and | are the right and left eigenstates 
corresponding to the complex eigenvalue of H% that has 
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the smallest positive real part among all eigenvalues in 
the non-rotating spectrum of H%. For a detailed review 
of this technique and related methods in non-hermitian 
Quantum Mechanics, see ref. [TT1I12) . The computational 
cost of this prescription scales exponentially with the 
number of particles. Since ng(r) — > hg{v) as |R| — » oo, 
and since there is a one-to-one correspondence between 
ng(r) and i>(re* e ) p~3l[T4], the complex energy of the LER 
Eg[ng] goes to E (not E), as |R| — > oo. Its lifetime C is 
given by (— 2Im(Eg))~ 1 , and for any finite |R|, 



Eg[ng] = £[n e ] - -C 1 [n g ] 



(3) 



where the resonance energy £ tends to E as |R| — > oo. 

To build a complex analog of Kohn-Sham DFT us- 
ing ng(r) as the basic variable, we first map the sys- 
tem of interacting electrons whose LER density is ng(r) 
to one of N particles moving independently in a com- 
plex "Kohn-Sham" potential uj(r) defined such that its 
N occupied complex orbitals {</>£ (r)} yield the interact- 
ing LER-density via n e {r) = T,tM' L \H^)\4' R ) ■ Thc 
complex Kohn-Sham equations are: 



hi — ei -h 2 - 2t- 
h 2 + 2t~ 1 hi - £i 



( Re(0! 



= , (4) 



where h\ = -icos(26»)V 2 + Re(uf(r)), and h 2 = 
isin(26>)V 2 + Im(uf(r)). The set of {ej and {tJ pro- 
vide the orbital resonance energies and lifetimes of the 
Kohn-Sham particles. 

Second, we write Eg[ng] as: 



Eg[n g ] = T e s [ng] + / dv n e (r)v{ve ie ) 



+E°[no] + E° 



[ng\ 



(5) 



in analogy to standard KS-DFT, and require: [tig] 
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T s [ng] and E„[ng] — e E H [ng], where T s [ng] and 
E H [ng) are the standard non-interacting kinetic energy 
and Hartree functionals evaluated at the complex densi- 
ties. Eq.pSlthen defines E^ c [ng]. The complex variational 
principle [12"] along with the assumption that the orbitals 
used to construct the density can be expanded in an or- 
thonormal basis leads to the Euler-Lagrange equation: 



5Eg[ng] 

Sng 



\i / drng(r) = 



(6) 



Performing the variarion in Eq. [5] and comparing with 
Eq.|4]leads to an expression for the Kohn-Sham potential 
that is again analogous to that of standard KS-DFT: 



J6\ 



e _ '%M(r) +i&.[ra e ](r) 



(7) 



where v* [n e ](r) = SE^ c [n e }/Sn e (r)\ LER . 

The simplest case where all essential aspects of this for- 
malism can be illustrated is a system of two interacting 



electrons moving in a one-dimensional potential such as 
the one depicted in the inset of Fig. [T] We study a Hamil- 
tonian where the electrons interact via a soft-Coulomb 
potential of strength A: 



using v(x) 



2 dx 



2 + v(Xi) 



A 



y/1 + (a?i - x 2 f 



2 

E 

3=1 



1 + e -2c(x+(-iyd) 



(8) 



Its 



parent potential v(x) = a(l - e- x2 ' b ) goes to a as x — > 
±oo, but v(x) goes down to zero at x ~ ±d. 

Exact solution via 2-electron wavefunction: The 

complex-scaled Hamiltonian Hg = H({xi} —> {xje }) 
was diagonalized with the Fourier Grid Hamiltonian 
(FGH) [15] and Finite Difference Methods. The numeri- 
cally exact ng(x) was calculated via Eq. [2] The complex 
density ng(x) depends on the value of 9 (see Fig. [I]), 
but for a large enough number of grid points the energy 
does not. In the complex-scaling method the resonance 
energies are precisely those that remain stationary as 9 
changes [T^]. Fig. [2] shows the energy for < A < 1. 

Exact KS solution : Two non-interacting electrons in 
the potential indicated by solid lines in Fig|3] have the 
same ng(x) as calculated above to one part in 10 6 (in 
the sense that the space integral of the square of the 
difference between their real or imaginary parts is less 
than 10 6 ). When ng(x) is set to integrate to the number 
of electrons (2, here), we verify this potential is given by: 



v!(x) = e-™^^-e H 
2y/ng{x) 



(9) 



where eh — ^Tr 1 ls the highest occupied complex orbital 
energy (in this case the only one), in exact analogy to 
real KS-potentials for bound 2-electron systems. With- 
out an explicit expression for E^ c [ng], however, the total 
energy cannot be calculated via Eq. [5] Related work by 
Ernzerhof P3] and physical intuition suggest that bound 
ground-state functionals are applicable here. They are, 
in any case, the most natural candidates. 

Exchange: Borrowing knowledge from bound 2- 
electron DFT, Eqs. [4] and [7j were solved employing 
Ei[ng\ = -\El\ng] = -\e-®E s [n e ]. The complex KS 
equations can be solved self-consistently with an initial 
guess for ng. Using the non-interacting complex density, 
the SCF calculations converged in 4-5 iterations. The re- 
sulting complex energies are plotted in Fig. [2] along with 
the exact results. For comparison, we also plot the re- 
sults from perturbation theory to first order in A. Thc 
two yield identical answers for the resonance energies, 
and extremely close for the lifetimes for all A in the range 
< A < 1. Thus, neglecting correlation, we find the av- 
erage error is ~ 14% for the real part and ~ 35% for the 
imaginary part of the total energy. We also compare with 
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FIG. 1. Different exact 2-electron complex densities when 
using different scaling angles. The model potential, v(x), used 
in this study is also shown (a = 4, b — 0.5, c = 4, and d = 2). 
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Re(E) 


hn(E) 


(N = 299) 0.27 


4.99895 


-0.0149586 


0.35 


4.99933 


-0.0144161 


0.43 


4.99962 


-0.0139792 


(N = 1299) 0.27 


5.00182 


-0.0161045 


0.35 


5.00198 


-0.0159848 


0.43 


5.00200 


-0.0159513 



TABLE I. Two-electron resonance energy values in the model 
Hamiltonian of Eq. [8] calculated via exchange-only DFRT. As 
the grid spacing decreases numerical dependence on 8 practi- 
cally disappears, (a — 4, b = 0.5, c = 4, d = 2 and A = 1) 

standard scattering calculations using the close-coupling 
equations under the bound state approximation [161 117) . 
The resonance energy is predicted by this method with 
an error of 22%, comparable to our DFRT exchange-only 
results. 

As in standard KS-DFT, total energies are given here 
by: 

N 

E e [n e ] = J2(e i -2irr 1 )+E e H3[ [ng] 

i=l 

- J drv e ux n e (v) (10) 

We point out that the ^-independence of the energy is 
preserved by the SCF procedure (see Table [I]). As the 
grid-size increases the dependence on 8 becomes negli- 
gible. This is important, because within a SCF DFRT 
calculation one is always solving the 1-body complex KS- 
equations. For these equations, one should be able to ef- 
ficiently use a large enough basis set or a fine enough grid 
to extinguish most of the numerical 9 dependence. Thus, 
this well-known drawback of the complex-scaling tech- 
nique [T8H2T)] is outdone by the benefit of never having 
to deal with iV-particle wavefunctions, but just 1-body 
(complex) densities. 

Correlation Potential : It is of interest to calculate the 
exact correlation potential, which we do by subtracting 




FIG. 2. The real and imaginary parts of Eg in the model 
Hamiltonian of Eq. [8] calculated exactly with complex scaling 
(thick solid), a first order correction to the non-interacting 
energy (dashed) , and our DFRT exchange-only self-consistent 
method (thin solid), (a = 4, b — 0.5, c = 4, and d = 2) 



the hartree-exchange contribution from the exact KS po- 
tential. The individual Hartree-exchange and correlation 
potentials are shown in Fig. [4] To interpret the features 
in these complex potentials it is useful to distinguish 
between two regions. As the interaction between elec- 
trons is turned on and A increases from to 1, the region 
around the central well is shifted up in the real part of 
the Kohn-Sham potential. This behavior is also seen in 
standard KS-DFT, and serves to shift up the position 
of the non-interacting orbital energies (in that case the 
real part of the orbital energies). However, both the real 
and imaginary part of the complex Kohn-Sham potential 
have a second region outside the central well that shows 
a dramatic oscillatory structure arising purely from the 
fact that the state is unbound. It is already known that 
the decaying oscillations in the tails of the complex LER 
wavefunction are governed by the lifetime of the reso- 
nance [H] . These oscillations serve to produce the correct 
assymptotic behavior in the interacting complex density 
and thereby give the correct interacting lifetime when 
this density is used in the functional. 

The analog of Koopmans' theorem does not hold in 
DFRT. Although the ionization energy of our 2-electron 
system is strictly zero, it is tempting to define Ig = 
Eq(N — 1) — E g (N — 2) and check whether it equals 
the highest occupied KS orbital energy. For the param- 
eters used in Figs.1-4, E g {N = 1) = 1.629 - 0.003i, 
Eg(N = 2) = 4.127-0.014i, but the exact KS eigenvalue 
is 2.065 — 0.006i. Clearly, DFRT provides an unambigu- 
ous prescription for the calculation of negative electron 
affinities. 

We are working on the implementation of DFRT to cal- 
culate the lifetime of molecular metastable anions. The 
method would also be applicable to molecules connected 
to metallic leads, as in molecular electronics. Ernzer- 
hof and co-workers have developed an approach for that 
purpose where complex absorbing potentials are added 
within a complex-DFT framework [il3j, [22] . However, we 
emphasize that the complex potentials in DFRT are the 
result of a variational calculation, and they are obtained 
self-consistently for the A^-electron system treated as iso- 
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FIG. 3. The real and imaginary part of the complex Kohn- 
Sham potential for the LER of 2 soft-Coulomb interacting 
electrons in the model potential, Eq[§] The dashed lines are 
the real and imaginary part of the complex-scaled parent po- 
tential v(x). (0 = 0.35, a = 4, b = 0.5, c = 4, d = 2, and 
A = l). 
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FIG. 4. The individual contributions to the Kohn-Sham po- 
tential from Hartree-exchange and correlation. (8 = 0.35, 
a = 4, b = 0.5, c = 4, and d = 2) 



lated, rather than added to the Hamiltonian from the 
start to model an open system. 

In addition, DFRT should be applicable to study shape 
and Feshbach resonances in low-energy electron scatter- 
ing processes of growing interest in biological 
systems [26"H25] . atmospheric sciences, lasers, and astro- 
physics 

In summary, DFRT provides an unambiguous prescrip- 
tion for calculating negative electron affinities based on 
a complex-scaled version of standard ground state DFT. 
This complex-scaled version has been cast in a way that 
is analogous in practice to KS-DFT. Results on a model 
system suggest that the same machinery that has been 
developed for KS-DFT yields accurate resonance energies 
and lifetimes in DFRT. It remains to be seen if common 
approximations to Exc [ n ] a re able to capture the impor- 
tant effects that determine properties of real transient 
anions. A more detailed study of the complex density 
function and various DFRT identities is forthcoming. 
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